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Abstract 

The metric for plane gravitational waves is quantized within the 
Hamiltonian framework, using a Dirac constraint quantization and 
the self-dual field variables proposed by Ashtekar. The z axis (direc- 
tion of travel of the waves) is taken to be the entire real line rather 
than the torus (manifold coordinatized by (z,t) is RxR rather than 
Si X R). Solutions to the constraints are proposed; they involve open- 
ended flux lines running along the entire z axis, rather than closed 
loops of flux. These solutions are annihilated by the constraints at 
all interior points of the z axis. At the two boundary points, the 
Gauss constraint does not annihilate the solutions, because of the 
presence of open-ended flux lines at the boundaries. This result is 
in sharp contrast to the situation in the general, 3+1 dimensional 
case without planar symmetry, where the Gauss constraints do not 
contribute at the boundaries because the Lagrange multipliers for 
the Gauss constraints fall to zero at spatial infinity. In the planar 
symmetry case, the Lagrange multiplier for rotations about the in- 
ternal Z axis survives at the boundary. The constraints are found 
to annihilate the solutions when classical matter terms are added to 
the Hamiltonian (so that flux lines are terminated on the matter). 

*e-mail address: nev@vm.temple.edu 
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The holonomy matrices used in the solutions are generated by the 
(2j+l) dimensional representation of SU(2), where j may be any 
spin, not necessarily j = 1/2. In this respect the solutions resemble 
the symmetric states, or spin network states recently constructed by 
Rovelli and Smolin in loop space. The Rovelli-Smolin area operator 
for areas in the xy plane is constructed and applied to the solutions, 
with puzzling results. 

PACS categories: 04.60, 04.30 
1 Introduction 

The connection-triad variables introduced by Ashtekar IQ have sim- 
plified the constraint equations of quantum gravity; further, these 
variables suggest that in the future we may be able to reformulate 
gravity in terms of non-local holonomies rather than local field op- 
erators [0, ^. However, the new variables are unfamiliar, and it is 
not always clear what they mean physically and geometrically. In 
particular, it is not clear what operators or structures correspond 
to gravity waves. Although the quantum constraint equations are 
much simpler in the new variables, and solutions to these equations 
have been found [@, ^, it is not clear whether any of these solutions 
contain gravitational radiation. 

This is the second of a series of three papers which search for oper- 
ator signatures for gravitational radiation by applying the Ashtekar 
formalism to the problem of plane gravitational waves. Paper I in 
the series |^ constructed classical constants of the motion for the 
plane wave case, using the more familiar geometrodynamics rather 
than Ashtekar connection dynamics. Paper II, the present paper, 
switches to connection dynamics, carries out a quantization of the 
plane wave metric, proposes solutions to the quantum constraints, 
and proves that the constraints annihilate these solutions except at 
boundary points. Paper III, which is in preparation, will propose 
operator signatures for gravitatonal radiation and apply those op- 
erators to the solutions constructed in the present paper. 

The phrase "plane gravitational wave" is used in this paper as 
a shorthand description for a specific class of metrics. "Plane" 



2 



means the metrics possess two commuting spacelike Killing vectors. 
I choose X and y coordinates so that these vectors are unit vectors 
pointing in the x and y directions. 

kW =9,;k(y) =9y. (1) 

The words "gravitational wave" may be defined precisely in a num- 
ber of ways IP] . Szekeres [0] suggests a definition which has consid- 
erable intuitive appeal: the system possesses an orthonormal tetrad 
consisting of k^^^ and k*^^) plus two null vectors which are hyper- 
surface orthogonal. The two null vectors may be interpreted as ray 
vectors for left-and right-moving radiation along the z axis. The two 
hypersurfaces are the corresponding wave fronts. 

It is possible to choose the coordinates (z,t) so that these wave 
fronts have simple equations ct±z = constant, but I follow the usual 
philosophy in the quantization literature and leave (z,t) arbitrary. 
When the (x,y) coordinates are gauge-fixed, the system becomes 
simpler than the full four-dimensional case. If the (z,t) coordinates 
are gauge-fixed as well, then the system becomes too simple and 
does not illuminate the full case. 

I quantize the theory in section 2 and derive wave functional 
solutions to the constraints in section 3. These solutions resemble 
Rovelli-Smolin T'^ operators, strings of triad operators separated by 
holonomies Q , except that the solutions fill the entire z axis (they 
are open, rather than closed flux loops) and each triad is integrated 
over z so as to guarantee invariance under z-diffeomorphisms. 

Also, the solutions utilize SU(2) generators which are (2j+l)x(2j+l) 
dimensional, rather than the usual 2x2 Pauli matrices. A 2x2 gen- 
erator can be viewed as acting upon a single flux line, while a 
(2j-|-l)x(2j + l) generator can be viewed as acting upon 2j flux lines, 
which have been totally symmetrized to give a total spin state of 
spin j. These new solutions may be closely related to the symmet- 
ric states, or "spin network" states recently constructed by Rovelli 



and Smolin in loop space ||I0|- I consider the introduction of the 
(2j-M)x(2j-(-l) generators as the most important technical innova- 
tion of the present paper. 

In quantizing the theory, I follow earlier work by Husain and 
Smolin M with one signiflgant exception: I incorporate a sugges- 



tion due to Teitelboim ITTI and rescale the scalar constraint. In 
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order to understand the advantages of this rescahng, consider the 
standard Dirac requirement that a wave functional in the physical 
subspace must be annihilated by an arbitrary linear combination of 
the constraints: 

= [NH + N^H, + NgHgItA 

= [(NE^z) (H/E|) + N^H, + NgHg] 

= [(N')(Hs)+N^H, + NgHg]^. (2) 

H and are the usual Asktekar scalar constraint and generator of 
z diffeomorphisms; Hg is the Gauss constraint, the generator of ro- 
tations around the internal Z axis. The planar symmetry has been 
used to gauge away the generators of x and y diffeomorphisms and 
generators of internal X and Y rotations. is the usual Ashtekar 
momentum variable, the densitized inverse triad ee^. The seemingly 
innocuous rescaling from H to Hg in equation (Q) has profound con- 
sequences for the algebra of constraints. Consider in particular the 
scalar-scalar commutator UTT 



[H(x), H(x')] = [Hi(x)ggy(x) + Hi(x')gg'j(x')]9j<5(x - x') + ■ ■ ■ 

- [H.(z)(E|(z))2 + H.(zO(E|(z'))^]9.5(z - z') + -(3) 

The • • ■ indicates terms proportional to the Gauss constraints, omit- 
ted for simplicity because they do not affect the argument. The first 
line of equation (|^) gives the commutator in the general, 3-1-1 di- 
mensional case; the second line specializes to planar symmetry. For 
a consistent quantization, the right-hand side of each line must an- 
nihilate ip, since the left-hand side does. In the general case, for 
most choices of operator ordering of the constraints H and Hi, the 
gg^j factor occurs to the right of the Hi and prevents the latter from 
annihilating ip. (It is possible to find orderings for which the gg'^ 
factor occurs to the left, but these orderings do not permit Hj to 
be interpreted as a generator of diffeomorphisms in the ith direction 
13[.) Now rescale as at equation (0). The rescaling modifies the 



commutator. 

[Hs(z), Hs(z')] = [H,(z) + R,{z')]dJ{z - z') + ■ ■ ■ . 

The unwanted gg^^ or E| factors have disappeared. As a result, in 
section 2, I will be able to pick a factor ordering which makes the 
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constraint algebra consistent and simultaneously permits to be 
interpreted as a generator of diffeomorpliisms. 

This is non-trivial progress; but it is not yet enough to guarantee 
a consistent quantization. The right-hand side of equation (0) could 
contain Schwinger terms, c-number terms which spoil consistency 
because they do not annihilate ip. If one had a complete set of so- 
lutions, presumably one could check for the presence of Schwinger 
terms by computing every matrix element of the commutator on the 
left in equation (|l]). Unfortunately, no such complete set is avail- 
able. However, Kuchaf [jl4[ has demonstrated that Schwinger terms 
are absent in a model which is very close to the present case, yet 
still solvable. The model is a free massless scalar field on a flat 1+1 
dimensional manifold. This model is made generally covariant by 
adding "embedding variables", plus scalar and vector constraints 
which guarantee that the new variables are not dynamical. Kuchaf 
rescales the constraints so that the algebra of constraints gives equa- 
tion (|I]) for the scalar-scalar commutator. Since he has a complete 
mode expansion, he is able to calculate Schwinger terms, and he 
finds no such terms in the constraint algebra. (The Lie algebra of 
the conformal group acquires Schwinger terms, in a manner familiar 
from string theory; but this conformal algebra is not the same as 
the Dirac algebra of constraints. Only Schwinger terms in the latter 
algebra destroy consistency.) 

Since Hs = H/E^, and H is polynomial, the new scalar constraint 
Hs is rational. I view this complication as a price which must be 
paid; but it is a small price, considering what one gets in return. 
The constraints close, which they should as a matter of principle; 
and the quantum theory possesses the diffeomorphism invariance 
characteristic of the classical theory |jT5l . 



Nevertheless, at some point a price must be paid: how does one 
define the inverse operator (E|)-i)? All but one of the terms in 
the Ashtekar scalar constraint H contain a factor of E|; therefore 
classically E| completely cancels out of all but one of the terms 
in Hs = H/E|. I will assume this cancellation occurs also in the 
quantum theory, so that Hg contains only one term with an E| in the 
denominator. The numerator of this particular term vanishes, when 
acting on the solutions constructed in this paper. Therefore, so long 
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as the action of (E|)~^ on ip does not give an infinity, (E|)~^ is not 
a problem in the present paper, and one can postpone a complete 
investigation of the inversion problem. In Appendix C, I determine 
enough of the action of (E|)~^ to show that this operator does not 
give an infinity, at least when acting on the solutions constructed in 
this paper. 

In three spatial dimensions it is usual to place the boundary 
surface at spatial infinity. Bringing the surface at infinity in to 
finite points is a major change, because at infinity the metric goes 
over to flat space, and flat space is a considerable simpliflcation. In 
the present case (effectively one dimensional because of the planar 
symmetry) the space does not become flat at z goes to infinity, and 
nothing is lost by considering an arbitrary location for the boundary 
surface. The "surface" in one dimension is of course just two points 
(the two endpoints of a segment of the z axis). These points will 
be taken to be a finite distance from the origin. The result that the 
space does not become flat as z goes to inflnity will be established 
in section 4; but for now, note that this result agrees with one's 
intuition from Newtonian gravity, where the potential in one spatial 
dimension due to a bounded source does not fall off, but grows as z 
at large z. 

Since the solutions are constructed from open flux tubes, there is 
internal "Gauss charge" exiting through the surface points. In 3+1 
dimensions, this exiting flux has no effect on dynamics, since the 
three Gauss charge operators Hj in the Hamiltonian are multiphed 
by Lagrange multipliers which fall to zero at spatial inflnity. The 
situation is different in the planar case. The high symmetry allows 
two of the Gauss constraints to be solved and eliminated from the 
theory. Prom the discussion of the behavior at spatial inflnity, given 
in sections 4 and 5, the Lagrange multiplier Nq for the surviving 
Gauss constraint Hz = Hq is finite at infinity. This means that the 
quantum Hamiltonian annihilates the solutions at finite z, but not 
at the boundary, because (the multiplier Nq is finite there and) the 
Hq term in the Hamiltonian is non-zero when acting on the exiting 
fiux. 

At this point one can think of (at least) four ways to proceed. (1) 
Ignore the difficulties at infinity and study the waves at finite z only. 
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(2) Insert classical matter sources proportional to 5{z — Zb), Zb the 
boundary points, into the Gauss constraint Hg, so that the modified 
constraint annihilates the solutions. The exiting fiux lines can be 
visualized as terminating on this matter. (3) Use second quantized 
matter sources : introduce Fermion fields into the fundamental La- 
grangian and terminate the fiux tubes on Fermion fields |TH]. (4) 
Search for closed fiux line solutions. 

The solutions (3)-(4) are the most satisfying in principle. Even in 
practice, presumably enough is known about Fermions to implement 
procedure (3) using the ansatz of section 3 for the gravitational 
part of the solution, although the calculation would be technically 
intricate. Procedure (4) probably requires a radically new ansatz; 
see the comments at the end of section 3. 

In this paper and the succeeding one I shall use procedures (1) 
or (2). These two procedures are identical in spirit: both ignore 
the difficulties at infinity, and study the waves at finite z. Section 
5 of the present paper constructs the matter sources needed for 
procedure (2). 

Procedures (l)-(2) have the virtue of simplicity, and indeed (1) 
is essentially the procedure employed in classical electromagnetism, 
where normally one studies a monochromatic plane wave. In elec- 
tromagnetism there are no problems with a Gauss charge, but the 
wave stretches from z = +oo to z = — oo and has infinite energy. 
The unphysical nature of this wave at infinity does not stop one 
from studying the wave at finite z and learning a great deal. 

In the electromagnetic case there is an analog of procedures (2)- 

(3) , the inclusion of matter. If one wishes, one can eliminate the in- 
finite energy by confining the radiation between two parallel metal 
plates. In practice, this procedure is never invoked until one wishes 
to discuss more advanced topics such as waveguides or Fresnel rela- 
tions at conducting boundaries. The introduction of metal bound- 
aries at any earlier point would be a distraction. I take the same 
point of view here: it is important to introduce realistic matter at 
some point; but initially the focus should remain on the radiation. 

The rest of this paper is organized as follows. Section 2 quan- 
tizes and regulates the theory. Section 3 proposes an ansatz for the 
solutions, and proves that the Hamiltonian annihilates the solutions 
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at all finite points away from boundaries. Behavior at boundaries 
is largely ignored in section 3, and in particular the surface terms 
in the Hamiltonian are ignored. Section 4 constructs these surface 
terms. Section 5 returns to the solutions of section 3 and shows 
that the Hamiltonian (now with surface terms included) does not 
annihilate the solutions at boundaries. Section 5 also implements 
procedure (2): matter sources are added at the boundaries so that 
the Hamiltonian does annihilate the solutions everywhere. Section 
6 discusses generalizations and directions for further research. 

Appendix A explains the sign conventions I have adopted in re- 
ducing the Ashtekar formalism from its covariant, four- dimensional 
form to 3+1 dimensions. The convention normally adopted in the 
literature has unintended consequences which I believe many 



physicists would prefer to avoid. 

My notation is typical of papers based upon the Hamiltonian 
approach with concomitant 3+1 splitup (except for the sign con- 
ventions explained in Appendix A). Upper case indices A, B, . . .,1, 
J, K, ... denote local Lorentz indices ("internal" SU(2) indices) 
ranging over X, Y, Z only. Lower case indices a, b, . . ., i, j, . . . 
are also three- dimensional and denote global coordinates on the 
three-manifold. Occasionally the formula will contain a field with a 
superscript (4), in which case the local Lorentz indices range over 
X, Y, Z, T and the global indices are similarly four-dimensional; or 
a (2), in which case the local indices range over X, Y (and global 
indices over x, y) only. The (2) and (4) are also used in conjunc- 
tion with determinants; e. g., g is the usual 3x3 spatial determinant, 
while '■^•'e denotes the determinant of the 2x2 X, Y subblock of the 
triad matrix e^. I use Levi- Civita symbols of various dimensions: 
ctxyz = exYZ = exY = +1- The Minkowski metric convention is 
'7TT = —1- The basic variables of the Ashtekar approach are an 
inverse densitized triad and a complex SU(2) connection A^. 

El = v^el; (4) 
[E1,A^] = -Mix-x')5l5l (5) 

The local Lorentz indices are vector rather than spinor. Strictly 
speaking the internal symmetry is 0(3) rather than SU(2), gauge- 
fixed to 0(2) rather than U(l). 
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2 Quantization 



This section begins with a brief discussion of the simphfications 
which ensue when the (x,y) coordinates are fixed, and follows this 
with a description of the Hamiltonian for plane waves. Then the sys- 
tem is quantized in a standard way, by replacing certain fields with 
functional derivatives. Issues such as factor ordering, regularization, 
and closure of constraints are discussed. 

The high degree of symmetry associated with the two Killing vec- 
tors allow Husain and Smolin to solve and eliminate four constraints 
(the X and y vector constraint and the X and Y Gauss constraint) 
and correspondingly eliminate four pairs of (E^, A^) components p[. 
The 3x3 matrix then assumes a block diagonal form, with one 
1x1 subblock occupied by E| plus one 2x2 subblock which contains 
all E^ with a = x,y and A = X,Y. The 3x3 matrix of connections 
A^ assumes a similar block diagonal form. None of the surviving 
fields depends on x or y. 

After these simplifications, the total Hamiltonian reduces to a 
linear combination of the three surviving constraints, 

Ht = N'[eMNEMEN(E|)~"^eABA^Ay + eMNEM^^b] 
+iN^E^Ffb 

-iNG[5,E|-eijEfAJ] + S.T. 
= N'Hs + N^H, + NgHg + S.T., (6) 

where 

K = d.A^ - 6nqA,^A^. (7) 
Hs,Hz, and Hq are the surviving scalar, vector, and Gauss con- 
straints. Strictly speaking these are Hamiltonian densities; for sim- 
plicity I have suppressed an integration over the z axis. S.T. denotes 
surface terms evaluated at the two endpoints on the z axis. The de- 
tailed form of these terms will not be needed until section 4. In 
equation (|^) the Lagrange multiplier N' and the scalar constraint 
Hg are rescaled versions of the usual Lagrange multiplier N and 
Ashtekar scalar constraint H, as at equation (^. 

In the general, three-dimensional case, often one chooses a wave 
functional depending on the A's, = ip[A]. Since the A's are con- 
nections, one is led to consider holonomies, and then a loop space 
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representation 0, In the planar case, only is a connection, 
and there is no special incentive to use any A except A^ in the 
wave functional. I shall choose a wave functional depending on this 
connection plus the fields 

% = [E|±i%]/y2, (8) 

where a = x,y. These fields are eigenstates of the surviving gauge 
invariance 0(2) or U(l) generated by Hq. Alternatively, one could 
use A^ eigenstates in the wave functional, but the weight one fields 
in equation (|^) will be signifigantly more convenient because they 
are densities and can be made diffeomorphism invariant simply by 
integrating them over z. 

The one connection field in the wave functional, A^, will be in- 
corporated into holonomy matrices in the usual way, 

M(z2,zi) =exp[i^'A,2Szdz], (9) 

M may be visualized as a flux line extending along the z axis from 
zi to Z2. There is an explicit factor of i because Sz is the usual 
(Hermitean) generator. 

Now The wave functional depends on the commuting set [Ej_,A^] 
which I will call the Q set, and does not depend on the conjugate 
P set [A;^, E^]. One may obtain a representation of the commuta- 
tors (^ in the usual way, by writing the P's as functional derivatives. 

At = h5/5E- 

= -h6/6Al (10) 

(If the pattern of ± signs seems strange, note that the two dimen- 
sional Kronecker delta in equation (^) has only off-diagonal elements 
when expressed in terms of 0(2) eigenstates: S±zp = +1.) 

Tentatively I choose to factor-order the Hamiltonian with P's to 
the right, Q's to the left, so that the constraint generates dif- 
feomorphisms along z. (More precisely, the constraint generates 
diffeomorphisms after a piece proportional to A^Hq is added [|12].) 
This ordering has been adopted already in equation (|^). 

The Hamiltonian must now be regulated, since it contains prod- 
ucts of functional derivatives evaluated at the same point. These 
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products can act on the wave functional and produce undefined 
products of delta functions evaluated at the same point. There are 
no products of functional derivatives in the Gauss and diffeomor- 
phism constraints, so that only the scalar constraint is in need of 
regularization. Every term in the original Ashtekar scalar constraint 
H contains a product of two functional derivatives, but the division 
by E| cancels one functional derivative from most terms, leaving 
behind a single functional derivative which no longer requires regu- 
larization. In fact the only term in Ht needing regularization is the 
(E|)~^ term in Hg. In the general, 3+1 dimensional case it is very 
difficult to find a regularization which does not violate diffeomor- 



phism invariance |jT8|, [T^, |l^ , but in the planar case one may employ 
a simple point-splitting regularization ||^. I arrange the three P's in 
the (Eg)"^ term along the z axis as follows. 

(E|(z))-ieABA^(z + 6)A,^(z - e) + {e-. - e) (11) 

As mentioned in the introduction, this term annihilates the solutions 
"0 constructed in the next section; consequently the exact form of 
this term matters little in the present paper, and a lengthy discussion 
here about the details of regularization would be pojntless. However, 
I do need to know enough about the action of the (E|(z))~^ operator 
in this term to prove that this operator is finite when acting on the 
solutions if). As shown in an appendix, the E^ operator can be 
inverted, and proved to be finite, provided the Sz SU(2) generator 
in the holonomy never has a zero eigenvalue. To see the reason for 
this restriction, note that equation (y) plus equation (|T0|) imply 

E|(z)e(zi+i, z, Zi)M(zi+i, Zi) = -ihSzMe(zi+i, z, Zi). (12) 

9 is a double Heaviside or square wave operator 

e(zi+i, z, Zi) = 6'(zi+i - z)6'(z - Zi). (13) 

If the Sz in equation ([12D is zero, then E^ has an identically zero eige- 
nalue and cannot be inverted. Therefore, I must assume Sz belongs 
to the (2j-|-l)x(2j+l) dimensional representation, withj half-integer; 
or, if j is integer, none of the Sz matrices in any holonomy takes on 
the value zero. 

There is one other context in which the (E|)^^ operator plays a 
role. Suppose I wish to prove that the constraint algebra closes, in 
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the following restricted sense: if I order the constraints on the right 
(in a typical constraint commutator such as equation (|I|)) so that 
P's are to the right, Q's are to the left in each constraint, and carry 
out the commutator, then the constraints on the left in equation (|I|) 
are also ordered correctly, with P's to the right, Q's to the left, and 
operators such as g'j absent or to the left. Although this is not full 
closure (nothing is said about possible Schwinger terms), restricted 
closure is a non-trivial result. Teitelboim's rescaling allows restricted 
closure. Since the proof of this statement is straightforward but 
lengthy, with the result presumably known already to Teitelboim, I 
will not go through the entire proof, but will mention only the step 
in the proof that is affected by the presence of the (E|)^^ operator. 
In order to commute this operator with operators also present 
in the constraints, I use 

[(El)-\A,^]=h5(z-z')(E|)-^ (14) 

This is the naive commutator one obtains by sandwiching the com- 
mutator of equation (|]) between two factors of (E|)^^. 

Although the double operator on the left in equation ( p!4D looks 
dangerous, it is well-defined on the solutions studied in this paper 
provided no Sz vanishes. This is enough to prove restricted closure 
valid when the constraints act on the solutions of this paper. 

The fact that the constraints close, in the restricted sense, with 
the operator ordering chosen, is an additional confirmation that my 
earlier, tentative choice of operator ordering is acceptable. I mention 
that the constraint algebra also closes for another popular choice of 
the P's and Q's: choose the P's to be the connections, and order 
them to the right. 

3 Solutions away from Boundaries 

This section proposes an ansatz for the solutions, then verifies that 
the Hamiltonian Ht annihilates these solutions at all finite points 
of the z axis away from boundaries. Strictly speaking it is not nec- 
essary to prove Ht'?/' = Ht(N', N^, Nq)'?/' = 0. (In fact after the 
surface terms are included, HtV' =(ADM energy) ^, which is non- 
zero.) Rather, one must prove Ht(^N', (5N^, (5Ng)'?/' = 0, where the 
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5N are arbitrary small changes in the Lagrange multipliers N. Small 
changes correspond to infinitesimal changes in coordinates, which 
should not make any physical difference. In this section, for sim- 
plicity of notation, I use N rather than SN. The distinction between 
N and SN becomes important at boundaries, where a given N can 
approach a constant, forcing the corresponding SN to approach zero. 

I return to the Hamiltonian, equation (|^) and break it up into 
eigenstates of 0(2) by writing out the components of the Levi-Civita 
tensor, 

e_+ = -e+_ = i, (15) 

while being careful to contract every + index with a - index, for 
example eMNAiviBN = e_+A+B„ + ■ ■ ■. When this is done one finds 
that half of the terms in H contain at least one factor of A^^ = 
hS/SE'^, hence half the terms would annihilate a wave funtional 
containing only E*^ fields. I try the ansatz 

V^(n;j) = n r"^'dzie(zi+i-Zi)M(z,+i,z„)E^^(z„)S+ X 

i=l "'^0 

xM(z„ z,_i) ■ ■ ■ E!_^ (zi)S+M(zi, zo). (16) 

Since this ansatz contains no Ei^ fields, half the terms in B.Tip vanish, 
including the term involving l/E^. Put aside the Hq term for the 
moment; it will require special treatment. The remaining terms all 
contain Fj^. 

Before evaluating the action of this operator on ip, note that each 
E(zi) is integrated over Zj (to preserve diffeomorphism invariance, as 
mentioned earlier), and the limits of integration have been taken at 
two finite but otherwise arbitrary points z = Zq and z = z^+i, rather 
than z = ±oo. As discussed in the introduction, nothing is lost by 
considering an arbitrary location for the surface. In fact something 
is gained: one no longer has to worry about the convergence of the 
dzi integrals in equation (p^. Because of these limits of integration. 



any theta function depending on zq or Zn+i in equation (|l^) may be 
taken to be unity. 

If the matrices S+ in ip were the usual 2x2 Pauli lowering opera- 
tors, then ^|J would vanish identically (except for the n = 1 case). I 
shall take S+ to be the (2j + l)x(2j + 1) representation, however, so 
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that i/j does not vanish identically, unless n > 2j. Similarly, the Sz 
in the holonomy, equation (y), is (2j + l)x(2j + 1). 

The surviving terms in (Ht — Hg)V' ^ill proportional to F^: 

(Ht-Hg)(z)V^ = -i(N' + N^)E^F+(z)^ 

oc &[d,{6tlj/6&) + iA^Sij/S&] 

= /" dzn ■ ■ ■ dzj^(zj+i - Zj)^(zj - Zj_i) ■ ■ ■ dziM ■ ■ ■ ME!^j(z) x 
j=i 



E---5(zj-z)E^-Hz)[9.j+iA,^|x 



j=i 

X [^(zj+i - Zj)^(zj - Zj_iM(zj+i, Zj)S_M(zj, Zj_i)] ■ ■ • + IBP. 



IBP denotes surface terms at zj = zq and Zn+i, which appear when 
we turn d^S into —d^jS and integrate by parts. All these surface 
terms are killed by the theta functions, except one at z^ = Zn+i 
and one at zi = zq. I put the IBP term aside for the moment 
and continue working on the main term. The d^^ acting on the M's 
produces a commutator — i[Sz,S4.]A^ which just cancels the iA^(z) 
term. This leaves terms in which d^j acts on the theta functions. 
These terms are 

(Ht-Hg)(z)V' = X:---E!:j+'(zj+i)M(zj+i,Zj)5(zj-z)E!:j(z)x 

j=l 

X - Zj)^(zj - Zj-i) + ^(zj+i - Zj)5(zj - Zj_i)] ■ ■ ■ + IBP, 

where delta functions with arguments depending on z^+i or zq are 
understood to be zero, and theta functions with these arguments 
are understood to be unity. Then the (j+l)st term in this sum 
contains a +5(zj_(_i — z)5(zj+i — zj) term which exactly cancels the 
— 5(zj — z)5(zj_,_i — Zj) coming from the jth term. (Ht — Hg)?/^ then 
collapses to the surface term IBP. 

(Ht-Hg)(z)V = ih(N' + N^)M(z„+i,z,+i)S+E-(z„+i)5(z-z,+i) X 

X J dZn_lM(Zn+l, Zn_l) ■■■ J dZi ■ ■ ■ M(zi, Zq) - 

-ih(N' + N^) J dz,M(zn+i,z„) ■ ■ .y dz2M(z2,zo) x 



14 



xS+E^i(zo)M(zo,zo)(5(z-zo). (17) 

For general choice of the E, (Ht — Hg)(z)'?/' is pure surface term, non- 
zero only at the boundaries z = zq and Zn+i- Note the delta functions 
imply (N'+N^) should be evaluated at the boundaries, and recall the 
remark at the beginning of this section that the Lagrange multipliers 
N are really small changes 6N. In the next section I shall require 
both N' and to approach constants at the boundaries; therefore 
the corresponding 6N vanish. These surface terms will turn out to 
be harmless. 

The calculation of Hq^ is very similar to the calculation just 
given. Cancellations occur at internal points, but terms survive at 
the boundaries. From equations and (p!OD, the 5zE| term in Hq 
produces Sz terms: 

- i9,E|M(zj, Zi) = h[SzM5(z - Zj) - MSz5(z - Zi)]. (18) 

As in the previous calculation, at points away from boundaries the 
factors of Sz commute with S+ and cancel the action of the other, 
eijEj term in Hq. At boundaries there is a surviving surface term: 

Hg(z)7/> = h[NG5(z - z„+i)Sz^ - Ng5(z - zo)^Sz]. (19) 

This term would vanish if Nq vanished at boundaries; but Nq does 
not vanish (or go to a constant) in the planar symmetry the 
next section will show. 

There is still the possibility that the Hqi/j term, equation (|19D , 



and the terms from Hx — Hq will be canceled by the surface terms 
in the Hamiltonian, which have the same support but have been 
omitted up to now. When surface terms Hgt are included and the 
Hamiltonian density is integrated over z, the constraint equation 
will have the following form. 

= [r^'dzHT(z)+Hst]V' 

= J dz[Hx(z) + (5(z — Zq) and S{z — Zn+i) terms]-?/'. (20) 

The next section will calculate H^t , and section 5 will check for 
cancellations. 
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It is possible to construct an additional set of solutions by re- 
placing every E!iS+ in equation ( [T6| ) by E!^S_. I have not been 
able to mix both El and E'^ in the same wave functional, however. 
The problem comes from a term proportional to N'[E'^Fj^ — E'i^F"^^] 
in Ht- The minus sign between the two F's destroys the delicate 
cancellations previously obtained between neighboring zj and zj+i 
terms in ijj. This lack of cancellation means that it is difficult to 
implement procedure (4) in the Introduction (procedure (4): search 
for closed flux solutions) using a wave functional of the present type 
(a string of E*^ and EJ^ fields separated by holonomies). Since a 
closed loop involves a trace, ip would have to contain equal numbers 
of E^ and E'^ fields (equal numbers of raising and lowering opera- 
tors S+ and S_). As suggested in the Introduction, implementing 
procedure (4) presumably requires a radically new ansatz. 

Is it possible to insert an arbitrary function f(zi ■ ■ ■ Zn) into ip7 In 
the linear case such a function (or its Fourier transform) is present 
in the wave functional and determines the spectral content of the 
wave packet. Here, however, one would obtain numerous dzi terms 
at the integration by parts step, and B.^ip would no longer vanish. 
Note that not even HsV' would vanish. This is signifigant because 
Husain and Smolin p| propose j = ^ f-dependent solutions which 
are annihilated by Hs (and Hq). They then transform to loop space 
in order to satisfy the remaining constraint. This trick does not 
work unless Hsip vanishes. 



4 Surface Terms 

The literature contains many discussions of surface terms. The clas- 
sical discussions based on metric dynamics [^, have been up- 
dated recently to systems with the most exotic boundary conditions 



22| . Within the framework of connection dynamics alone there are 
at least three recent discussions |T^, . Nevertheless, the pla- 

nar case has enough twists and turns to make another discussion 
interesting, as well as necessary. 

The typical discussion of surface terms begins with two assump- 
tions. The first assumption (algorithmic assumption) justifies the 
need for the surface terms and provides an algorithm for calculating 
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these terms. Since the algorithm typically requires detailed knowl- 
edge of the behavior of the basic fields on the boundary, a second 
assumption (asymptotic assumption) is needed to supply this be- 
havior. 

I utilize an algorithmic assumption proposed by Regge and Teit- 
elboim When the Hamiltonian (more precisely, the classical 

Lagrangian) is varied to obtain the equations of motion, 5[d^ (field)] 
variations occur which must be replaced by 5 (field) variations. The 

may be removed by integration by parts, which in turn gives rise 
to a surface term containing ^(field). The total variation is required 
to vanish (surface term as well as volume integral). The vanishing 
of the volume integral gives the classical equations of motion, of 
course, but the surface term does not vanish, in general, unless one 
adds a compensating surface term to the original Hamiltonian. 

To see how the algorithm works in practice, let us apply this 
method to the present Hamiltonian Ht, equation d). Vary Ht and 
look for 5 [dz (field)] terms. Every field strength F conttibutes a 
term of the form ES^JA. Also, the Gauss constraint Hq contributes 
a dzSE^ term. After integrating by parts to remove from each 
5 (field), 

S J dzHx = volume integral + [N'ewNE^^A^ + iN"E^(5A^ 

-^5E^A^] ■ (21) 

I have used a result from the 3+1 decomposition, that the Lagrange 
multiplier Ng is also the four-dimensional connection field 



P!7| . One can try to cancel these surface terms by adding to the 
original Hamiltonian J^^'^^ Ht a surface term Hgt of the form 

J dzHT + Hst = / dzHT - [N'eMNE^A^ + iN^E^A^^] |^;'+^ . (22) 

At first glance this surface term does not seem to work. Its varia- 
tion does indeed give 6 A terms which exactly cancel the correspond- 
ing 6 A terms in equation (0); but 6B.st seems to be missing an Af 
term needed to cancel the corresponding Af term in equation (pT]); 
and (5Hst contains 5E^ , 6N^, and 6N' terms which are unwanted. I 
could supply the missing A^ term by adding to Hst an additional 
term +iE|A^ ; but again, this term generates an unwanted 5A'^ 
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variation. For now I will eschew terms. (We shall see later that 
they are unnecessary; equation ( p2D is correct as it stands.) 

At this point I need the asymptotic assumption to show that the 
"unwanted" SE^ terms have exactly the right asymptotic behavior 
to cancel the term in equation (^iD , while the quantities N' and 

approach constant values which are not subject to variation at 
the boundaries (hence 6N^ and 6N' vanish there). To formulate an 
asymptotic assumption, one falls back on classical experience and 
adopts for the quantum case the simplest asymptotic behavior that 
works in the classical case for the system of interest. In the case 
of a three-dimensional system with bounded sources, for instance, 
the asymptotic assumption is flat space at infinity. In the one- 
dimensional radiative case, the analog of a bounded source is a wave 
packet or packets, located inside the region zq < z < Zn+i, and zero 
at the boundaries (zq, Zn+i). From Newtonian gravity it is too much 
to hope that the metric for this system is flat at the boundaries, but 
Szekeres ^ uses hypersurface orthogonality of the null tetrads plus 
the field equations to show that an assumption of conformal flatness 
in the variables (z,t) is always possible . To make this idea precise, 
I introduce Szekeres' parameterization for the plane gravitational 
wave metric. 

ds^ = e^[dx^e^ cosh W + dy^e~^ cosh W - 2dxdy sinh W] 

+ e°-^/2{[-(N')' + (N")']dt2 + 2NMzdt + dz'}. (23) 
Conformal flatness at the boundaries means 

N'^1;N"^0, (24) 

so that the (z,t) sector of the metric assumes a e°^^/^[— dt^-|-dz^] = 
— 2e°~^''^dudv form, flat except for a scale factor. Equation (p^ ) is 
already a considerable simplification, since it means N' does not 
have to be varied in Hst, and all surface terms may be dropped. 

It should be possible to simplify the metric still further at the 
boundary, since gravitational wave degrees of freedom vanish there. 
In the linearized limit, the fields B and W defined in equation (p^) 
are amplitudes for the two polarizations of the gravitational wave. 
This result suggests that in the exact theory B and W should be 
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taken as wave packets which vanish at the boundary. 

B^O;W-^0. (25) 

Co nsistent with assumption (^3]), the exact classical equations of 
motion for B and W are hyperbolic. (This result is straightforward 
to prove but requires a digression into the classical equations of mo- 
tion and formulas from paper I, and I relegate the proof to Appendix 
B.) 

Now only the fields exp(A) and D remain at the boundaries. In 
the linearized limit, D = and exp(A) = 1, while in quadratic order 
exp(A) obeys a parabolic, rather than hyperbolic equation with the 
linearized energy as source 

2(expA)„,+[(B,)2 + (W,)2 + (ttb)' + (vrw)'] = 0. 

TTx is the momentum conjugate to X. Thus exp(A) is a good candi- 
date for long-range scalar potential, but not for wave packet behav- 
ior. D is gauge-sensitive and could be anything, from the arguments 
given so far . 

However, D turns out to vanish at the boundaries. The argument 
was given in paper I, but is perhaps worth repeating here, because 
it does not use the usual asymptotic assumption. Paper I is based 
on geometrodynamics, so that there is no B.q constraint, and the 
Hs and constraints are expressed in terms of the fields A, B, D, 
W and their conjugate momenta. To determine D at the boundary, 
I need only the part of Ht which is independent of B, W, and N^, 
therefore does not vanish at the boundary: 

6[ I dzHT + Hst] = 6 1 N'[2(exp A)„, -(expA)„D„ + ■ ■ •] + 5R,t 

= volume terms + [2(5expA),z — ((5exp A)D,z 

-(expA)„5D + 5H3t] l^f'' (26) 

where ■ ■ ■ denotes terms which vanish at the boundary or contain no 
derivatives. Obviously the -|-2(5exp A),^ term can be canceled by 
inserting a — 2(exp A),^ term into Hst; but what term can I insert to 
cancel the 5expA and 6D variations? It turns out there is no such 
term, and I must take 

D^O (27) 
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at boundaries, to eliminate variations that cannot be canceled. 

To prove that there is no such term, one may try to construct such 
a term until frustration sets in; or one may construct a formal proof 
as follows. Suppose such a term exists; call it f. The functional f = 
f[expA,D] occurs in Hgt, and its variation cancels the D-dependent 
IBP terms coming from the volume term: 

= -(5expA)D„-(expA)„ 5D + j dz'{5e^(z')[5f(z)/5e^(z')] 

+5D(z')[5f(z)/5D(z')]. 

This implies 

5f(z)/5e^(z') = 5(z-z')D,; 
(5f(z)/5D(z') = 5{z-z'){e\. 

This is impossible because the two second functional derivatives are 
not equal, QED: 

5H{z)/5e^{z')5T){z") = 5{z-z")dJ{z-z'). 

The asymptotic behavior of all fields except D was determined 
by the usual asymptotic assumption (which is: the quantum fields 
have the same asymptotic behavior as the classical fields). D, on the 
other hand, was determined solely by the algorithmic assumption, 
without invoking the asymptotic assumption at all! No cancelling 
Hgt term can be found; therefore the Lagrangian formulation is not 
consistent unless D is set equal to zero. 

So far I have used the algorithmic and asymptotic assumptions 
to find the asymptotic behavior of the usual metric fields (more pre- 
cisely the fields B, W, ■ ■ ■ used by Szekeres to parameterize the met- 
ric). It is not a good idea to apply the algorithmic and asymptotic 
assumptions directly to the Ashtekar triad and connection fields, 
because their classical behavior is poorly understood. However, it 
is now straightforward to obtain the asymptotic behavior of the 
tetrads, triads, and connections, since these fields may be expressed 
in terms of the metric fields. 

1 start with the tetrads. In the (x,y) sector the general formulas 
for the tetrads as functions of B, W, • • • is quite comphcated, but 



20 



the formulas are needed only at the boundaries, where the metric is 
diagonal and the expressions simplify considerably. 

e±^5^exp(A/2±i0) (28) 

0, the angle of rotation around the Z axis, is the new field which 
appears on switching variables from the 3-metric (four independent 
metric fields) to the triads (five independent triad fields). The clas- 
sical behavior of this field is 

0(z, t)^constant, (29) 

so that 50 = at boundaries. The expression for Hgt is a scalar 
under internal Z rotations, hence cannot depend on 0. Without loss 
of generality I may take 0^0, therefore. 

The (z,t) sector of the tetrad is simple enough that one can 
write down expressions true for all z, then take their limit as z 
— ^-boundaries: 



e^ 



= N'yg^^^exp(-A/4); (30) 

el = yg^^-.exp(-A/4); (31) 

= N^v^^,-0; (32) 

ej = 0. (33) 

(These tetrads embody the standard choice of gauge for Lorentz 
boosts along the Z axis; for the advantages of this choice, equa- 
tion (H), see Peldan 0.) 

Given these tetrads, one can now determine the E fields from the 
relation 

E\ = ee[. (34) 

The result is 

E| = expA; (35) 
E^ ^ 5^exp(A/4). (36) 

Relation (|35|) is correct for all z. 

Now one may compute the connections A[ by computing the four- 
dimensional Lorentz connection uj^^ , then expressing the A[ in terms 
of the a;^''. 

^" = [^i[ja] + ^[ai] - f^a[ij]]e''ej^ (37) 
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where 

a[ja] =eiM[ajef -aaef]/2. (38) 

(Since the fields u and are four- dimensional, they should carry 
a superscript (4) which I suppress for simplicity.) Equation (|37|) 
simplifies considerably at boundaries, or wherever the tetrad matrix 
becomes diagonal. 



ul'^[e,K{d^ef-d-4)]^'^\ (39) 
Using these equations, one finds 

- 0; (40) 
cuj^ ^ exp(A/4)9texp(-A/4); (41) 
cuf ^ -5Cexp(-A/4)a,(expA)/2; (42) 
a;™ ^ 5^xp(-A/4)at(expA)/2; (43) 

- 0; (44) 
uoj'^ exp(A/4)(9,exp(-A/4), (45) 

where B, C = ± only. These are the only Lorentz connection com- 
ponents needed to calculate the A[. The latter follow from the 
equations given in Appendix A: 

-iexp(A/4)9texp(-A/4); (46) 

PAB _ ^ , ,ZC • TB 
LrA^ - -CBZC^^a " l^a 

^ -eBaexp(-A/4)9,(expA)/2 

-i(5fexp(-A/4)9t(expA)/2; (47) 

^ -iexp(A/4)(9,exp(-A/4). (48) 

Note the last equation: in contrast to the situation in three dimen- 
sions, the Lagrange multiplier for the Gauss constraint does not go 
to zero at the boundaries. 

Since the asymptotic behavior of the triads and connections de- 
pends on the asymptotic behavior of the exp A field, I take a moment 
to investigate the asymptotic behavior of the latter. It is relatively 
easy to extract this behavior from the geometrodynamical equations 
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derived in paper I. After a canonical transformation to the A,D,B,W 
parameterization, the Lagrangian assumes the form 

L = TTAAt + TToDt H Ht 

= 7rAA,t +7rDD,t H 

-N'[2(e^)„, -(e^)„ D„ -e'WTrB + ■■■], (49) 

where ttx is the momentum conjugate to X, • • • denotes terms de- 
pending on B, W, and (which do not contribute near boundaries) 
as well as the Hst term (which may be ignored since its variation does 
not give the classical equations of motion). Setting the scalar con- 
straint equal to zero gives a parabolic equation for e^ which can be 
solved provided the e~^7rA7rD term is known. From the variation of 
the Lagrangian with respect to ttd, 

D,t +N'e-^7rA = 0. 

Hence from equation ( P7[ ) for D, tta vanishes at the boundary. When 
this result is inserted into the scalar constraint, (expA),^^ is found 
to vanish; consequently expA at boundaries behaves like the New- 
tonian scalar potential in one dimension. 

exp A^ciz + C2, (50) 

where ci and C2 are time-dependent functions. Moreover, 2ci is the 
ADM energy [Q, therefore does not vanish, in general. 

I now have the asymptotic behavior of all fields which occur in the 
surface terms, and can return to the question raised at equation (pOl) : 
do the surface terms from Ht, equation ^T]) , cancel the surface terms 
from the variation of Hst, equation (p2D ? The total surface term 
obtained on variation is 

S J dzHT + 5Hst = volume term + [-i^E.^A,^] \^^;' -[N'eMN^E^A^] 

-> volume term + [-i5(e'^)(i)(A/4),z ] 

-[eMN5(5Mexp(A/4))(-l)eNbexp(-A/4)9,(expA)/2] 
= volume term + [e'^5AA,^ /4] - [e'^^AA,^ /4] 
= volume term. (51) 

On the first line the first square bracket comes from Ht and the 
second from Hst. N^^O terms have been dropped; and the 5A^ 
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terms from the two brackets have canceled out. On the second hne, 
the Kronecker part of A^, equation (^Tj), has been dropped; it 
cancels in a later sum over indices. Note in the first bracket the 

Gauss term not only survives to infinity; it is actually needed at 
infinity, to cancel the (5E^ term from the second bracket. From the 
last line of equation (^) all surface terms have canceled from the 
Euler-Lagrange variation, QED. 

I have also experimented with a more complex Hst which includes 
an A^ term: 

H^, = aiHst + a2iE|A,^ 

where ai and a.2 are adjustable constants. When varied, the new 
surface term generates not only e'^^AA,^ terms, of the same form 
as equation (^Tj), but also new terms of the form 6{e^),z- Requir- 
ing both types of variation to vanish independently places two con- 
straints on the two constants aj, and the only solution is ai = 1, a2 = 
0. There are no Af terms in Hgt 



5 Solutions at Boundaries; Classical Matter 

The surface terms Hst derived in section 4, equation (P2|), can now be 



used to answer the question raised in section 3: does Hst cancel the 
integration by parts boundary terms obtained when Ht is applied 
to ipl The expression which must vanish is 

( J dzHx + Hst)V'. (52) 

More precisely, must be invariant under infinitesimal transforma- 
tions generated by the Hamiltonian. If Ht = Ht[N', N^, Nq], and 
similarly for Hst, then the expression which must vanish is 

( J dzHT[5N', 5N^ 5Ng] + Hst[0, 0, 5Ng]^, (53) 

and it must vanish for arbitrary choice of the small quantities 6N. 
Note the zero arguments in equation (|53|): two of the N's are fixed at 
the boundary by equations (p^), and only Nq = Af can be varied at 
the boundary. However, from equation (|22D there are no Nq or Af 



terms in Hst. The surface term in expression (Rl) is identically zero. 



The integration by parts surface terms from Ht remain uncanceled. 
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The introductory section suggested several ways to modify the 
Hamiltonian or wave functional so as to eliminate the surface terms 
and obtain a solution to the constraints. The remainder of this 
section will implement procedure (2) suggested in the Introduction: 
add approximate Fermionic matter terms to the Hamiltonian and 
include the Fermion fields in the wavefunctional, so that the flux 
lines terminate on Fermion fields and the modified wave functional 
is annihilated by the Gauss constraint. In order to describe quan- 
titatively how much flux exits at each boundary, I introduce some 
notation. Let m^ be the Sz eigenvalue of the homotopy matrix to the 
left of E'i(zi)S+ in ■0. When previously suppressed matrix indices 
are exhibited, 

= [M(z„+i,z„)E-(z,)S+---E-_Hzi)S+M(zi,zo)]'"^,„ 
= M-"^^E-(S+)'X-i ■ • ■ E^i(S+)"^UM-?,„. (54) 

Evidently the matter at z^+i must have Sz — mn, while the matter 
at zq must have Sz = mo. 

The added Fermionic matter terms may be either quantized (op- 
erator fields) or classical (c-numbcr functions). Suppose first the 
matter terms are quantized; take them to be (say) the Hamilto- 
nian for a free spin 1/2 Weyl Fermion. "Free" means that covariant 
derivatives are replaced by ordinary derivatives and spatial tetrads 
are replaced by unity, but Lagrange multipliers are left unchanged, 
so that the matter terms appear as free-field additions to the usual 
scalar, vector, and Gauss constraints. Then the algebra of con- 
straints will continue to hold. (Note that the structure constants 
of this algebra are independent of the gravity fields, because of the 
rescaling suggested by Teitelboim; hence the interacting Fermion 
and free Fermion terms obey the same algebra.) One can make 
ip into a Gauss scalar by sandwiching ip between initial and final 
Fermion wavef unctions. These are constructed by multiplying to- 
gether enough spin 1/2 Fermion fields to get the required total Sz 
of mo or mii_|_i. (At this point each Fermion in the product should 
be given a distinguishing "color" index, to avoid difficulties with 
Fermi statistics.) If the spin 1/2 Fermi coordinates are used to 
construct these wavef unctions, rather than the canonical Fermi mo- 
menta, then these initial and final wavefunctions are not densities 
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and do not have to be integrated over zq and Zn+i. This construc- 
tion is straightforward enough, but for consistency I would have to 
muhiply the gravitational wave functional by a wave functional for 
the Fermi field (or rather by a product of such wave functionals, one 
for each color). This seems unnecessarily elaborate. 

I therefore turn to classical matter. Here there is no hope of 
getting a closed algebra of constraints, but closure of the matter 
terms is not required for consistent quantization anyway. Consider 
for example the scalar-scalar commutator, which produces a linear 
combination of the vector and Gauss constraints. Schematically, 

[Hs.Hs] = H, + Hg. 

Even if the scalar constraints on the left have c-number matter 
terms, they certainly will not survive commutation and appear on 
the right. However, recall that these commutators are smeared with 
small changes 6N'i and 6N'2- These changes vanish at the bound- 
aries z = Zb, hence if I take the c-number corrections to be propor- 
tional to 5(z — Zb), they will disappear from the commutator algebra. 
(The [Hg,Hg] commutator does not give trouble, even though the 
smearing function 6Af is finite at boundaries, since this commuta- 
tor vanishes identically.) I therefore propose the following c-number 
modification of the Gauss constraint. 

Hg^Hg + hmo(5(z - zo) - h.mj{z - z^+i). (55) 

There are analogous, delta-function modifications of the Hs and H^ 
conmstraints, but their explicit form will not be needed. 

I must now return to sections 2 through 4 of this paper, and check 
that all boundary conditions and gauge choices adopted in those 
sections continue to hold in the presence of the matter. In section 
2 ("Quantization") I eliminated four constraints by adopting the 
gauge due to Husain and Smolin p. I must check that this gauge 
choice still eliminates these constraints (the Hx, Hy diffeomorphism 
constraints and the Hx, Hy Gauss constraints) when Fermionic mat- 
ter is present. Although I do not quantize the Fermion matter, for 
orientation it will be helpful to have at hand the Hamiltonian for 
the free spin 1/2 Weyl Fermion. 

Hp = -N'ujOyr] - N^^ud-^T] + iAluja^r]/2. (56) 
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T]^ is a spin 1/2 Weyl Fermion and uj-q is its conjugate momentum. 
In this section I reserve Roman letters A,B, ■ ■ ■ from the beginning 
of the alphabet for spin indices (A,B,- ■ ■ = ± 1/2), and continue 
to use Roman capitals I, J,- ■ ■ from the middle of the alphabet for 
the Lorentz indices (I, J,- ■ ■ = X,Y, or Z). Thus is shorthand for 
the Pauli matrix I note that the diffeomorphism constraints 

Ha with a = (x,y) in equation (|56D contain (9a. Since this vanishes 
because of the symmetry generated by the spacelike Killing vectors, 
there are no matter corrections to Hx or Hy in equation (|56D , and I 
will assume that there are likewise no c-number matter corrections to 
Hx or Hy in my model. In equation ( ^61) the Hx and Hy constraints 
contain I = (X,Y). If mo is positive (say), I can build the initial 
Fermion wavefunction entirely out of spin 1/2 Fermions having spin 
+1/2, while imposing the constraint that the spin -1/2 rj and uj 
fields vanish. Then the terms for I = (X,Y) are constrained to 
vanish and there are no matter corrections to Hx or Hy. Taking 
my cue from this example, I assume that the c-number Fermions 
used to construct the initial wavefunction in my model are entirely 
polarized in the mo direction, and similarly for the Fermions in the 
final wavefunction, so that there are no c-number matter corrections 
to Hx or Hy. 

Husain and Smolin fix the x and y diffeomorphism gauges and 
the X and Y Gauss constraints by choosing 

E| ^ ^ 0, (57) 

where a = x,y, I = X,Y, and ~ as usual denotes a gauge condition 
which is to be imposed only after carrying out all commutators. Now 
if one commutes the Hamiltonian (including matter terms) with the 
constraints of equation (|^, one finds that these constraints are 
conserved in time only if further constraints hold: 

kl ^ Ai^O; (58) 
k\ ^ N^^O. (59) 

These are exactly the additional constraints obtained by Husain and 
Smolin, as one would expect, since the added matter terms are c- 
number and do not change the equations of motion for the tetrad 
and connection fields. Further commutation of the constraints (|58|) 
with the Hamiltonian yields no new constraint. 
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At this point I can treat the constraints in equation ( |57D and 
equation (|58| ) as four pairs of second-class constraints, and ehminate 
them using the Dirac bracket procedure. Since the four pairs are 
just canonical coordinate-momentum pairs, the "Dirac brackets" for 
the remaining E and A fields are just the usual Poisson brackets 
for these fields. After setting the constraints in equation ( ^7D and 
equation (^) strongly equal to zero, I reproduce the Smolin-Husain 
result that the gravitational part of the Ha and constraints are 
strongly equal to zero, for a = x,y and I = X,Y. Since the Fermionic 
part of these constraints also vanishes, all four of these constraints 
vanish strongly, and I may set = A]: = for a = (x,y) and I = 
(X,Y): I no longer need the equations obtained by varying N'^ and A{ 
in the Lagrangian. At this point I have eliminated four constraints 
and four coordinate-momentum pairs, the same number as Husain 
and Smolin. Their gauge choice continues to work in the presence 
of c-number matter terms. 

In section 3 (" Solutions" ) it is easily verified that the modified Hq 
of equation ( |55| ) now annihilates the solution ip of section 3: com- 
pare equation ( |55D to equation (]T9|). In section 4 ("Surface Terms") I 
must review Szekeres' argument ||^ that the metric may be brought 
to conformally fiat form at infinity. Szekeres uses the Newman- 
Penrose spin coefficient formalism and shows that a certain spin 
coefficient a must vanish. From this it follows that the metric de- 
composes into two disconnected 2x2 subblocks, a (z,t) subblock and 
an (x,y) subblock. Then he uses the existence of two hypersurface- 
orthogonal null vectors to bring the (z,t) subblock (globally) to 
a conformally fiat form. The Newman- Penrose a happens to be 
a combination of the connections which appear in equation (p8|) . 
Therefore the Szekeres argument remains valid, because the Husain- 
Smolin gauge remains valid. 

Since the added matter terms are being treated as c-number 
"background" matter, they are not to be varied when obtaining the 
equations of motion; consequently, they contribute no integration 
by parts surface terms; no matter terms need be added to Hgt. The 
argument that D— >-0 at the boundary therefore continues to be valid. 
Only the argument that exp(A) is linear at the boundaries, equa- 
tion (0), must be modified, since the scalar constraint now contains 
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matter terms. All other results in section 4 continue to hold. This 
completes the check that the boundary conditions and gauge choices 
of sections 2-4 continue to hold in the presence of c-number matter. 

I comment briefly on the effect of including interactions in the 
matter terms. For example, adopt the quantum operator Hamilto- 
nian of equation with ordinary derivatives replaced by covariant 
derivatives and with tetrad fields restored. Because of the covari- 
ant derivatives, the Fermion Hamiltonian now contains connec- 
tion terms. If one tries to impose the gauge choice of Husain and 
Smolin, equation {^7\), the secondary constraints of equation (|58|) 
are replaced by constraints of the form ^ (bilinear combination 
of Fermion fields). This is a manifestation of the fact that in a the- 
ory with minimally coupled Fermions and no "curvature squared" 
terms, the gravitational connection acquires a torsion, with the tor- 
sion non-dynamical and bilinear in the Fermion fields. Because of 
the torsion, the Newman- Penrose a no longer vanishes. Evidently 
in a theory with interacting Fermions, even the kinematics (gauge 
choices) would have to be rethought from the beginning. 

6 Directions for Future Research 

The Rovelli- Smolin symmetric states |T^, or spin- network states 
r[7] are functionals of loops. The loops 7 are flux lines of spin 
j which join and separate at trivalent vertices; the vertices are 
Clebsch-Gordan coefficients. Symmetric states and the states '?/'[E, A] 
introduced in the present paper obviously resemble each other, but 
the exact relation between the two is not straightforward to estab- 
lish. In three dimensions a functional Fourier transform connects 
r[7] to '?/'[E, A]. This transform was established by demanding that 
certain T" operators defined in (A, E) space have simple analogs in 
loop space 0. When three dimensions are replaced by one, the T° 
algebra changes dramatically, since most of the A^ are no longer 
connections, and the Fourier transform must be rethought from the 
beginning. This is beyond the scope of the present paper. 

However, the Rovelli- Smolin construction does suggest the follow- 
ing generalization of the solutions of section 3. In ip the (2j + l)x(2j+l) 
S+ raising operator multiplying each E^^ is proportional to a Clabsch- 
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Gordan coefficient 



(SK)m'm = -VJ(J + 1) < Jm' I IKjm > . (60) 

(Since the coefficient on the right obeys Condon-Shortley conven- 
tions p6|, the Sk on the left must be normahzed hke a Condon- 



Shortley vector operator, which means S±i = =f(Sx ± iSY)/V^- 
One can generalize equation (|60|) to 

(SK)m'm (const.) < jW | IKjm >, (61) 

where j' may be j ± 1 as well as j. Simultaneously, the (2j + l) 
dimensional holonomy M to the left of S+E^^: must be replaced by 
a (2j' -|- 1) dimensional M. Note the only property of the Sk matrix 
used in deriving the solutions of section 3 is 

[Sz,Sk] = KSk. (62) 

Written in terms of Clebsch-Gordan coefficients, this is an obvi- 
ous identity which generalizes immediately from equation (^) to 
equation (|6lD, so that the matrices of equation (pT| ) also generate 
solutions. 

Rovelli and Smolin have constructed a gauge-invariant, diffeomorphism- 
invariant, and regulated area operator. I have applied this operator 
to the solutions of the present paper, with puzzling results. This 
work is described in Appendix D. 

Husain and Smolin have constructed a large number of solutions 
which are more general than the solutions presented here, in the 
sense that the wave functional involves both S+ and S_ [§. These 
solutions use only j = 1/2 and satisfy the constraints because of the 
special properties of the Pauli matrices. Some (very) preliminary 
work indicates that these solutions cannot be generalized to the 
case j > 1/2. The j = 1/2 sector may be the best place to look 
for a lowest energy state, however, and the Husain-Smolin solutions 
deserve further study. 



A 3+1 Phases 

When the four-dimensional covariant formalism (with real Lorentz 
connection '-^•'a;^'') is rewritten as a 3+1 dimensional canonical for- 
malism (with complex connection A|), one must make three choices 
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of phase. Two of these choices are straightforward and require httle 
discussion: choose the phase of the Levi-Civita tensor by choosing 
the sign of exxYz; and choose between self-dual or anti-self dual 
complex connection by choosing a phase 6 = ±1. 

2G(^)AiJ =(^) + i6{e'\j2e^^y,ruf\ (63) 

The duality relation is 

i5(6"KL/2eTXYz)(^)Af =(^) Ai^ (64) 

The explicit factor of exxYZ simplifies later formulas and guarantees 
that S is independent of choice of phase for the Levi-Civita tensor. 

The third choice of phase occurs when the connection ^^^A^^ with 
two space-time Lorentz indices is replaced by the connection A[ with 
a single spatial Lorentz index. For Lorentz indices M,N = space- 
space, the equation connecting the one- and two-index connections 
is 

(')Ar = a6MNsA^/2, (65) 

where a = ±1 is the arbitrary choice of phase. For Lorentz indices 
T,M = time-space, the equation connecting one-and two-index con- 
nections follows from equation (BBl) plus duality, equation 



(^)A™ = -ia6Af/2. (66) 
The standard proof that the Ashtekar formulation of gravity is 



equivalent to the usual formulation |23, |13 starts from a gravita- 
tional lagrangian L expressed in terms of the real connection ^'^^cu; 
then a pure imaginary topological term is introduced, so that L can 
be rewritten as a function of the complex connection ^'^-'A; then equa- 
tions (|65D -(|66D are used to eliminate *^^^A and bring the Lagrangian 
to pq — Ht canonical form. The dependence of the final result on 
the above phase conventions exxYZ, and a is as follows: 

87rL= - ia(5E^A^,t+NaeMNsESE?^FL/2 

+ ia^N'^EJ^jFl + ia^E^DnAf . (67) 

The phase exxYZ has completely dropped out of the final result. 
One can read off the canonical momentum from the coefficient of 
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A^,t. Since the momentum depends on a6, so does the equal-time 
commutator after quantization: 

[A^,EN] = -a5h<C5(z-z'). (68) 



F and Dn in equation (pTf ) are the three-dimensional field strength 
and covariant derivative respectively. 

Fmn = <9mA^ — (9nA^ — aesMNAj^A^G; 
D.A^ = 9„Af - aeNMRA^^^A^^G. (69) 

So far there is zero reason to prefer any particular phase for 
ctxyz; and little or no reason to prefer any particular phase for 6. (I 
have chosen +1 for both these phases.) What about a? The usual 



choice in the literature appears to be +1 |T^. However, note that 
traditional conventions for the matrix elements of the 3x3 J = 1 
SU(2) generators Ss imply 

Aa[Ss]MN = +eMSNA^, 

with S in the middle rather than at the end, which suggests a = —1 
in equation (p5|). If a = -|-1 is actually the unconventional choice for 
this phase, one should get unconventional results in equations which 
are non-linear in the fields A^. Indeed, consider equations (|5P[). 
With a = +1, these are not the traditional definitions of field 
strength and covariant derivative in SU(2) gauge theories. In this 
appendix I have left a and 6 as arbitrary phases in the final formu- 
las, so that authors may make their own choice, but in the body of 
the paper I have chosen e-rxYZ = 6 = +1, a = —1 . 



B Equations of Motion for the B and W Fields 

This appendix verifies the assertion made at equation (p5D that the 
classical equations of motion for the B and W fields are hyperbolic. 
From paper I, the Lagrangian written in terms of Szekeres fields B, 
W, D, A is 

L = TTB cosh WB,t+7rwW,t-N'e^[(B„ cosh W)2 + (W„)V2 

-N'e-^[(7rB)' + (vrw)']/2 - N^Ittb cosh WB„ +7rwW„ ] + (70,) 
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where ■ ■ ■ indicates irrelevant terms depending on A and D and their 
conjugate momenta only. I will prove the assertion for the B field. 
The proof for the W field proceeds similarly and is slightly easier, 
because the equations contain no cosh W factors. 
The variation of L with respect to ttb and B gives 

= cosh WB,t -N'e-'^TTB -N^ cosh WB„; 
= -(7rBCOshW),t+[N'e^B„(coshW)2]„ 

+N^[7rBCOshW]„. (71) 

One can obtain a wave equation for B in the usual way, by solving 
the first of these equations for ttb and inserting the solution into the 
second equation. 

= -[e^(l/N')(B,t -N^B„ )(coshW)2],t +N'e^B„, (cosh W)^ 

+ N^[e^(l/N')(B,t -N^B„ )(cosh W)2]„ + ■ ■ ■ , (72) 

where ■ ■ ■ are terms containing no second derivatives of B. Note 
that to establish the hyperbolic character of wave equation (|72[), 
one needs only the second derivative terms ||2^. Carring out the 
differentiations in equation ([7^), keeping only the second derivatives, 
one gets 

= -e^(l/N')(B,tt-N^B„t)(coshW)' + N'e^B„,(coshW)' 
+N^e^(l/N')(B,t, -N^B„, )(cosh W)2 + ■ ■ ■ 
= e^(cosh W)2(l/N'){-B,tt +2WB,^, +[(N')' - (N^)']B„, } + ■ ■ ■ 
= aB,tt +bB,2t +cB,22 H . (73) 



The condition for a hyperbolic equation is b — 4ac > [^] , which 
translates to [e^(cosh W)^]^ > 0, QED. Note that it is necessary to 
divide through by N' and e"^ at several points in the above proof. If 
these quantities vanish, the metric of equation (^3]) becomes degen- 
erate, since the determinants of the (z,t) or (x,y) subblocks vanish. 



C Inversion of E| 

This appendix computes the action of 1/Ez on enough functionals 
to show that 1/E| is bounded when acting on the solutions ip of 
section 3. 
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To invert an operator, such as the momentum operator p in quan- 
tum mechanics, one might try the following (naive) procedure. Ex- 
pand an arbitrary state in a complete set of momentum eigenstates. 
On each eigenstate the action of p gives a constant, 

pexp(ikx) = hkexp(ikx). (74) 

Then define the inverse as one over the constant: 



[1/p] exp(ikx) = [1/hk] exp(ikx). (75) 

This definition has the advantage that p[l/p] gives unity, as it 
should. The problem is that k can vanish. 

When the spectrum of the operator is continuous (as is the case 
for p) one can try defining 1/k at k = 0, say by replacing k— >-k ± ie. 
There is no point in discussing this option here, since the spectrum 
of E| appears to be discrete: from equation ( plSj) its action upon a 
holonomy gives a factor of Sz- Therefore vanishing eigenvalues are a 
problem, and the conclusion given in section 2 appears unavoidable: 
holonomies with Sz = must be excluded. Once the kernel of E| is 
excluded, E| can be inverted using the method sketched above for 
p: find the eigenfunctions, and define the inverse as one over the 
eigenvalue. 

Strictly speaking there is not one operator E^ (z) to invert, but 
an infinite number, one for each value of z. The eigenvectors also 
depend on z, if the operators depend on z: O(z)0(z) = A(z)0(z). 
The holonomy times a B function is a z-dependent eigenvector: 

E|(z)[e(zi+i,z,Zi)M(zi+i,Zi)] = -ihSz[e(zi+i,z,Zi)M] (76) 

9 is the square wave function defined at equation (plsl). The holon- 



omy M is defined at equation (^, and I have used equation (plOl). 
The inverse of equation ([76|) is 

[l/E|(z)][e(zi+i,z,Zi)M(zi+i,Zi)] = i(hSz)-Me(zi+i,z,zOM] (77) 

Next consider the action of l/E^ on the solutions ip of section 3. 
These solutions have the form 



= M(z,+i, z„)X(zn) ■ ■ ■ M(z2, zi)X(zi)M(zi, zq). (78) 



34 



The X are operators in the Lie algebra of SU(2) and are independent 
of . Zn+i and zq are the boundary points, and for simphcity I 
suppress the integrations over the remaining dzi. Multiply ip by the 
following partition of unity. 

1 = [e(Zn+l, Z, Zn) + 6(Zn, Z, Z^-l) H h 0(z2, Z, Zi) + e(zi, Z, Zq)]. 

(79) 

Then 

[1/E|]^ = [l/E|][e(z,+i,Z,Zn) + --- + e(z2,Z,Zi) + e(zi,Z,Zo)]^ 
= i(hSz)"^6(Zn+l, Z, Zn)M(Zn+l, Zn)X(Zn) ■ ■ ■ M(zi, Zq) 
+ ■ ■ ■ + M(Zn+l, Zn)X(Zn) ■ ■ ■ i(hSz)-'e(zi, Z, Zo)M(zi, 

where equation ([77|) was used on the second line. Provided no Sz 
vanishes, this action is completely finite, QED. 



D The Area Operator 

The classical expression for the area of a surface S with normal 
is 

A = f dVv'EfE^'nanb. (81) 

s 

In order to regulate this product of operators, Rovelli and Smolin 



29| evaluate the E operators at two slightly separated points a{s) and a(t) 
on a small loop a, then insert holonomies between the separated E 
operators to maintain gauge invariance. In order to obtain diffeo- 
morphism invariance, they square A and doubly- integrate over the 
surface. The result is 

A'= f dV / dV|(l/2)c2T^N,nb|, (82) 

S s 

where 

T'^^[a](s,t) = Tr[M(s,t)Ef(a(t))((TVic)M(t,s)E!^(a(s))((TVic)], 

(83) 

and the holonomies M are integrated along the small loop a. Rovelli 
and Smolin renormalize all their Pauli matrices with (an arbitrary 
real factor c times) a factor i, as in equation (|83|). This is done pri- 
marily to remove an explicit factor of i from the holonomy; however. 
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Rovelli and Smolin then have to contend with two factors of i in the 
definition of the area operator. They have removed these from equa- 
tion (P^ ) by using absolute value bars rather than a factor (-1). As 
a result, the a^O limit of equation (^) contains jEjEj*! rather than 
simply Ej E^. The absolute value bars are signifigant and I shall 
return to them later. 

First, however, one must ask how the above, 3+1 dimensional 
construction must be modified to fit the planar case. It is possible 
to define area operators for areas lying in the xz, yz, and xy planes. 
Since the ip are functionals of only a single connection A^, however, 
only the xy area operator oc E^ is a simple functional derivative. 
The other two area operators multiply ip hy a, function, and there 
is no hope of a simple action. I will therefore concentrate on the 
xy area operator. The construction of this operator is immediate: 
replace the general holonomies in equation (^) by holonomies along 
the z-axis; replace both Ef by E^; replace the factor ic by my nor- 
malization, a factor of 2. 

Since the factors of i have now disappeared from T'^'^, there is no 
need for absolute value bars (or a factor (-1)) in equation (R3). The 

resultant A operator has a small a limit oc E| , the limit Rovelli 
and Smolin would have obtained had they used a factor (-1) rather 
than absolute value bars. 

Now for simplicity, and to avoid the complexities of the Penrose 
calculus, let A^ act upon the holonomy constructed from the j = 
1/2 representation. This holonomy is infinitely functionally differ- 
entiable, so that the point-splitting in equation ( ^3l) is overkill, and 
one may replace the holonomies in T'^'^ by unit matrices. Then the 
area operator collapses to 

A^ = [ydxdyE|]2. (84) 

Each E| brings down a factor of — ihS^ (or hr^ if one uses the Rovelli- 
Smolin normalization, with the factor of i hidden in cx cr^/i). 
Thus A^ has a negative eigenvalue: 

A^M = (-ihS^)2M 

= -{hS^fU. (85) 
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The eigenvalue (S^)^ is 1/4, only 1/3 of the Rovelli-Smolin value 
j(j+l) = 3/4, but the 1/3 can be accounted for by the shift from 
three spatial dimensions to one. The factor of (-1) is the important 
feature: it comes from the factor of i in the holonomy, and has 
nothing to do with dimension. This factor would be removed by the 
absolute value bars, but it is not clear that the absolute value bars 
are demanded by any fundamental principle of the theory. 

This difficulty can be made more striking. Note that in the planar 
case it is possible to take the square root of the operator and still 
get a diffeomorphism invariant result. The square root 



is also gauge invariant, since the only surviving gauge rotation (around 
the Z axis) leaves the Z index invariant. This operator has an imag- 
inary eigenvalue when acting on the holonomy. 
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